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ABSTRACT 

Recent near-IR (NIR) and X-ray observations of Sagittarius A*'s spectrum 
have yielded several strong constraints on the transient energization mechanism, 
justifying a re-examination of the stochastic acceleration model proposed pre- 
viously for these events. We here demonstrate that the new results are fully 
consistent with the acceleration of electrons via the transit-time damping pro- 
cess. But more importantly, these new NIR and X-ray flares now can constrain 
the source size, the gas density, the magnetic field, and the wave energy density in 
the turbulent plasma. Future simultaneous multi-wavelength observations with 
good spectral information will, in addition, allow us to study their temporal evo- 
lution, which will eventually lead to an accurate determination of the behavior 
of the plasma just minutes prior to its absorption by the black hole. 

Subject headings: acceleration of particles — black hole physics — Galaxy: center 
— plasmas — turbulence 



1. Introduction 

On 27 October 2000, the Chandra X-ray observatory detected a highly variable X-ray 
flare coincident with the position of Sagittarius A* (Baganoff et al. 2001). The transient 
lasted a couple of hours, with a peak luminosity ~ 45 times greater than the quiescent 
emission. Even more surprising was the realization that during the event, the X-ray output 
dropped abruptly by a factor of five in under ten minutes, recovering just as quickly. Light 
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travel-time arguments therefore place the source of this unusual radiation within a region no 
bigger than about 10 light-minutes or ~ 10 13 cm across. 

Sagittarius A*, a compact radio source at the Galactic center, is thought to be the 
radiative manifestation of a ~ 3.4 x 10 6 M supermassive black hole with the Schwarzschild 
radius r s = 2GM/c 2 ~ 10 12 cm, where G, M and c are the gravitational constant, the black 
hole mass and the speed of light, respectively. Earlier theoretical modeling of its spectrum 
and polarization properties (Melia, Liu, and Coker 2000, 2001; see also Falcke and Markoff 
2000 for an alternative model in which this emission is produced within a jet) had already 
anticipated an emission region within the inner ten Schwarzschild radii of a hot, magnetized 
Keplerian flow. 

Not long after Chandra's first detection of the X-ray flare, XMM-Newton followed with 
its own measurements, including the discovery of two unusually strong bursts a couple of 
years later (Goldwurm et al. 2003; Porquet et al. 2003). None of the previous X-ray 
satellites had the sensitivity and spatial resolution to identify such low-luminosity events at 
the distance of the Galactic center. Chandra and XMM-Newton now detect them at a rate 
of about one per day, most of which are usually weak and last tens of minutes. The best-fit 
photon index during the majority of these bursts is ~ 1.3, representing a flattening of about 
1 compared to Sagittarius A*'s spectrum in the quiescent state, which includes a significant 
contribution from thermal emission at large radii (Melia 1992; Baganoff et al. 2003). In 
addition, the 27 October 2000 event appears to show a soft-hard-soft spectral evolution. 

But the most intriguing X-ray flare of all may actually be the most recently detected 
event, in which an unambiguous modulation with an average period of 21.4 minutes was 
seen over the course of its ~ 3 hour duration (Belanger et al. 2005). The separation 
between the flux minimums actually decreases from about 25 down to 17.5 minutes as the 
flare evolves, corresponding to the passage of an emitting plasma in a Keplerian motion 
from a radius r ~ 2.9 rs to 2.4 rs- This region therefore appears to lie somewhat below 
the marginally stable orbit (MSO) for a non-spinning black hole; it would, however, be 
outside the MSO for a Kerr black hole with a large spin. The monotonic decrease of the 
X-ray period is strongly reminiscent of what was seen in near-IR (NIR) flares detected just 
a few years earlier, where an average period of 17 minutes was associated with a similar 
chirping behavior with the period decreasing from 23 to 13 minutes (Genzel et al. 2003). 
These two sets of observations — one ground-based in the NIR, the other at X-ray energies 
from space — support the view that we are probably witnessing the evolution of an event 
moving inwards through the last portion of the accretion disk inside or very near the MSO. 
The inferred radial velocity v r ~ 10 8 cm s _1 is consistent with an accretion driven by the 
magnetic viscosity of the turbulent Keplerian flow. 
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Not surprisingly, previous speculation on the underlying mechanism for these events (Liu 
and Melia 2002) focused on the view that such X-ray flares might be driven by an accretion 
instability. While this may still be true in light of all the more recent observations, the process 
by which the actual emission occurs is uncertain. However, the fact that the NIR spectrum 
(Eisenhauer et al. 2005) is much steeper than that of its X-ray counterpart excludes a direct 
extrapolation of the spectrum (see Liu and Melia 2001 and Fig. 2 below). The currently 
favored scenario is one in which the mm/sub-mm to NIR portion of the spectrum is due 
to synchrotron, whereas the X-rays are produced by synchrotron-self-Compton (SSC). (This 
constraint is empirically motivated, and is independent of whether the emission occurs within 
a disk or a jet.) It turns out that producing the right blend of physical conditions to fit both 
the NIR and X-ray flare emission (under the assumption that the two occur more or less 
simultaneously — a concept that is yet to be confirmed compellingly) is not trivial. In the 
next section we describe the observational constraints on the model parameters under this 
scenario. Our work with stochastic acceleration (SA) in producing Sagittarius A*'s quiescent 
spectrum (Liu, Petrosian, and Melia 2004, LPM04 hereafter; Liu, Melia, and Petrosian 2005) 
motivates us to consider a picture in which the flare itself is produced by a magnetic event, 
possibly driven by an accretion instability. A toy model of this acceleration and its fit to the 
flare spectra are described in § 3. § 4 summarizes the main results and discusses the model 
limitations. 



2. Observational Constraints on the SSC Model 

The NIR and X-ray flares in Sagittarius A* have peak luminosities as high as ~ 10 35 
ergs s _1 (Baganoff et al. 2001; Ghez et al. 2004). Several relatively long duration flares, two 
in the NIR (Genzel et al. 2003) and another in X-rays (Belanger et al. 2005), also displayed 
quasi-periodic modulations with a period decreasing from ~ 25 to 13 minutes as the flare 
evolved. Assuming Keplerian motion, this corresponds to a transition in radius from ~ 2.9rs 
to 1.9r£. Two flares have also been observed simultaneously in the NIR and X-ray bands, 
with a peak X-ray luminosity, respectively, 3 and 18 times the quiescent level of ~ 2 x 10 33 
ergs s" 1 and their spectroscopy still being processed (Eckart et al. 2004; Baganoff et al. 
2005). The first NIR spectroscopy was completed in July 2004: a power- law fit to the power 
spectrum vF v oc v~ a (with u the emission frequency) yields a = 2.2 ± 0.3 during the peak 
of the flare observed July 15, 2004, and a = 3.7 ± 0.9 during the rising and decay phases. 
For the flare of July 17, 2004, a = 3.5 ± 0.4 (Eisenhauer et al. 2005). If the NIR emission is 
produced via the synchrotron process, the radiating electron distribution N^) oc 7 _p (with 
7 the Lorentz factor) must have an index p > 7.4, suggesting an exponential cutoff of the 
electron distribution presumably dictated by the acceleration process, at 7 cr ~ (^nir/^b) 1 ^ 2 ) 
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where ub = eB±/m e c is the electron gyre-frequency, e, m e and B± are the electron charge, 
mass and the perpendicular magnetic field, respectively. The X-ray flares, on the other 
hand, often display a very hard spectrum with a ~ —0.7 ± 0.5 (Baganoff et al. 2001; 
Goldwurm et al. 2003). In the SSC (or in general IC) scenario, this requires a flat electron 
spectrum (p ~ 1.6). (This is different from a power-law commonly assumed, or a broken 
power-law distribution caused by radiative cooling.) At the longer sub-mm wavelength the 
flares usually have a much smaller flux increase above their quiescent value than the NIR 
flares, which suggests a photon spectral index a < 2.2, requiring a flattening of the electron 
distribution at lower energies. As we shall show below a fairly flat power-law spectrum with 
an exponential cutoff N(j) = A" 7~ p exp(— 7/7 cr ) can reproduce these observed spectra and 
is a natural consequence of a simple SA model. Because most of the observed NIR and 
X-ray emissions are produced by electrons near 7 cr , and for p < 0.4 the radiation spectrum 
is almost identical, we set p = —2, corresponding to a relativistic Maxwellian distribution, in 
what follows. The observed emission characteristics can then set strict limits on the model 
parameters, such as the source size R, the magnetic field B, the gas density n, and 7 cr . 

Before describing the acceleration model we discuss how the existing observations limit 
the possible ranges of these parameters. As mentioned above the very steep NIR spectrum 
requires the cutoff frequency of the synchrotron emission v CT = 3e-B_i_7 2 r /47rm e c = 1.7 x 
10 12 (7 cr /100) 2 ( J B ± /40 G) Hz in the sub-mm to NIR range (Eisenhauer et al. 2005). This sets 
an upper limit on 7 2 r -B_i_. The spectrum of the SSC photons is also very sensitive to B± and 
7 cr . In the SA model described below, the scattering rate is much higher than the acceleration 
and energy loss rates. The electron distribution is isotropic. We therefore consider the pitch 
angle averaged results. The solid and long-dashed lines in Figure 1 (left) represent three 
spectra produced by electrons with a Maxwellian distribution with the spectral indexes 
a = — din {yF v )/d\n.(v) at v — 1.4 x 10 14 Hz (q;nir) and v = 10 18 Hz («x) indicated by 
the label on the lines. The change in a as one moves about the B-^ cr plane is due solely to 
the dependence of v cr on B and 7 cr . For flares with a soft NIR spectrum (onir < 0) and a 
hard X-ray spectrum («x > 0), one can exclude the upper-right and lower- left portions of 
the B-^cr plane. 

The acceleration time r ac , which is independent of energy and equal to one quarter of the 
energy loss time at 7 cr in the SA model discussed below, must be shorter than or comparable 
to the flare rise time r r j S < 60 minutes, except that the source is strongly Doppler-boosted 
toward the observers. Because the luminosity of X-ray flares is usually smaller than that 
of the sub-mm to NIR flares, we can set the energy loss time equal to the synchrotron 
time T syn (7) = 9mlc 5 /Ae A B 2 ^ = r /7, where we have assumed that the source is optically 
thin. (The SSC cooling due to a radiation field with an energy density U syQ can be readily 
incorporated by replacing B 2 by an effective field B 2 S = B 2 + 8nU syn whenever there exists 
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an observational justification.) The dashed lines in Figure 1 (left) give three acceleration 
times. The left-hand side of the B-j cr plane can be excluded because r ac > 60 mins there. 
So to produce a flare with q;nir > 0.0, ax < 0.0, and r ris < 60 mins via synchrotron and 
SSC, the maximum Lorentz factor and magnetic field must be located within the central 
region of Figure 1. In principle, simultaneous NIR and X-ray spectroscopy during the flares 
can thus directly fix B and 7 cr . 

Let us now consider how simultaneous NIR and X-ray observations with good spectral 
information for both can in addition provide us with a measurement of R and the density 
n of the radiating electrons. The synchrotron luminosity due to an isotropic relativistic 
Maxwellian electron population can be estimated as follows (Pacholczyk 1970) 

ioe KfBW 

where M = 2N / y^ r is the total number of the accelerated electrons. For an isotropic radiation 
field, the corresponding X-ray emission produced via SSC is then given by 

r U syn 87r£ syn 

= 52xlo35 (iow) 2 (i^)"(l)" ergss "' (2) 



where Ub = B 2 /8tt and A is the magnetic field energy density and the surface area of the 
source, respectively, and £ syn ~ U syn cA. From these we can obtain A and N ' . For a uniform 
spherical source, A = AttR 2 and Af = AnR 3 n/3, we have 

R ~ 0,64 (l0 36 ergs s" 1 ) (l0 35 ergs s" 1 ) (loe) Ts ' (3) 

- - 4.6 x 10 6 ( J syn T ) " ( / SSC T ) ^ ( -*-) (^) " cm- 3 . (4) 
V 10 36 ergs s" 1 J V 10 er g s s J \A0G J \100J w 

This source size and density are typical values expected in the accretion torus of Sagittarius 
A* (LPM04), and, in general, depend on the source structure and geometry. The right panel 
of Figure 1 shows how these quantities may be read directly from simultaneous NIR and 
X-ray spectroscopic observations: From the NIR (solid) and X-ray (long-dashed) spectral 
indexes, one can locate the flare in this parameter space; The bottom axis then gives the 
magnetic field; With the NIR and X-ray flux density measurements, R can be read from the 
top axis; The left and right axes then give the density n of the accelerated electrons and the 
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turbulence to magnetic field energy density ratio / tur b, which we shall discuss in the next 
section. 

Flares are likely triggered by some plasma instability, or by changes in the dynamics 
of the accretion flow — for example, the dissipation of angular momentum is different above 
and below the MSO, which may lead to a strong gravitational dissipation and acceleration 
of electrons. It is important to note that the total number of energetic electrons required to 
produce these bright flares can set a limit on the mass accretion rate M and accretion time 
Taccr : T accr M > Afm p , where m p is the proton mass. If r accr ~ r ris , then 

*>„. 9xl „»(£)(^)- gl! -., (5, 

which is consistent with the value inferred from Sagittarius A*'s quiescent state spectrum 
and can account for the observed flares for a radiation efficiency of ~ 10%. 



3. Toy Model of Stochastic Electron Acceleration 

Let us now turn our attention to the SA of electrons and examine how this process meets 
the challenges posed by the above constraints. Electrons can be accelerated by turbulent 
plasma waves via the transit-time damping process and cyclotron resonances. The former 
dominates in plasmas where the gas pressure is higher than the magnetic field pressure (Yan 
& Lazarian 2004), and the pitch angle scattering rate is higher than the acceleration rate 
when the Alfven velocity va = B(A^nm p )~ 1 / 2 < c (Schlickeiser & Miller 1998), resulting in 
an isotropic electron distribution. The corresponding acceleration rate is proportional to the 
electron energy, giving rise to an energy independent acceleration time: 

T ac = 2 C 2 R f — , (6) 

which is determined by the turbulence generation length ~ R, f tUT b, an d (Miller, LaRosa, 
& Moore 1996 1 ). The SA can therefore be addressed by solving the following kinetic equation 
for the spatial and pitch angle integrated electron distribution function N (Petrosian and 
Liu 2004) 



dN _ d 
dt c?7 



d'-f N ( 47 2 t ; 



4 7 _ ] N 



N 



T, 



+ Q, (7) 



CSC 



1 Notc r ac = jm e c 2 /A, where the acceleration rate A (not to be confused with the source surface area 
defined above) is given by Equation 2.4a in the paper. 
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where the source term Q and N/T esc give the rate of particles entering and escaping from the 
acceleration site, which we identify as the emission region. As we shall see below, the energy 
density of energetic electrons is at least one order of magnitude higher than the magnetic 
field energy density. We therefore envision that the plasma is excited by certain internal 
instability. In such a scenario, the escape time T esc can be much longer than the duration of 
the flare. 

For a constant injection rate 0(7) — Qo$(l — Tin) with 7 in < 7 cr , where 

Tp _ 97T 2 mgcX/turb _ Qn (R U y 1 ( ' ft^b\ (R , 

7cr "4r ac " lQe^RB 2 " \r s ) \W cm^J \ 0.1 J ' [) 

the steady-state spectrum cuts off exponentially above 7 cr but is different from the Maxwellian 
one (Park & Petrosian 1995; Schlickeiser 1984): 

/°° /9 2r \ ^ 2 

x 5 - 1 (l+x) 3+<5 exp(- 7 x/ 7cr .)dx with 5=(- + j^\ -1.5 

(9) 

The left panel of Figure 2 shows this distribution for 7 in = 0.137 cr and for several values of 5 
as indicated. The solid line gives the Maxwellian distribution. For large T csc or small 5, most 
of the observed emissions are produced by electrons near or above 7 cr , consequently the slope 
of the spectrum below 7 cr is unimportant, and the Maxwellian distribution provides a good 
approximation of the accelerated particle spectrum. For iV( 7 ) = iVo7 _p exp(— 7 / 7cr ), our 
calculations show that nearly identical photon spectra can be produced for different values 
of p by adjusting N . 2 

There are therefore four model parameters: B, R, n, and /t ur b (or 7cr)- 111 the previous 
section, we have shown that the observed emission characteristics of NIR and X-ray flares 
set strict limits on these quantities. A detailed spectral fitting for individual flares can be 
used to measure them. The thick solid line in the right panel of Figure 2 fits the peak 
spectrum of the NIR flare of July 15, 2004, using R = 0.4 r s , B = 40 Gauss, n = 1.4 x 10 7 
cm -3 , and / tur b = 0.14 (or 7 cr = 75). Here the NIR spectrum alone determines £?7 2 r and 
Af, and the X-ray spectrum, which we assume to be hard and a factor of ~ 20 above 
the flux density of the quiescent state, constrains R and B^ r . This model is designated 
by the black dot in Figure 1. If the NIR flares always have very soft spectra, say with 
a > 1.4, then the magnetic field associated with these flares will fall into the narrow range 
16(r syn /60 mins)- 2 / 3 G < B < 120 G, and 7cr will be between 53 and 140(r ac /60 mins) 1 / 3 . 



2 Note that should the escape time be short, i.e. when 8 ~ (2T ac /T csc ) 1 / 2 is large, there could be a huge 
low energy electron population for 7i n <C j C r- Equation (8) shows that /turb — 1 is required to accelerate 
electrons to GeV energies. 
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The source is self-absorbed in the sub-mm and longer wavelength range. For the source 
size inferred from the NIR flares the expected sub-mm flux will be ~ 0.3 Jy, which would be 
difficult to detect decisively. However, SA accompanying these flares in a large volume can 
produce energetic electrons, which produce strong mm and longer wavelengths emission in 
situ or in the process of escaping toward larger radii (Zhao et al. 2004). It is also interesting 
to note that the acceleration time r ac is 27 mins for this model, and the energy flux associated 
with the cascading Kolmogorov turbulence is ~ fll^R 2 vpJJ B ~ I0 35 ergs s" 1 (Miller et al. 
1996). These values are fully consistent with the flare observations. 

In the right panel of Figure 2, the thin solid line corresponds to a model with B = 38 
Gauss, which fits the rising and decay phase spectrum of the NIR flare of July 15, 2004. 
The dashed line fit corresponds to another model with B = 90, which fits the peak flux of 
the X-ray flare of 27 October, 2000. All other model parameters have the same values as 
those described above. These results suggest that the variations in flare characteristics may 
be attributed to changes in the magnetic field. However, we emphasize that the optically 
thin NIR emission only depends on the total number of energetic electrons, while the X-ray 
emission also depends on the source area. To demonstrate this effect, the dotted line shows 
a spectrum produced with the area of the emission region increased by a factor of 4 (i.e., 
R = 0.8rs- Note M is fixed). As expected, the NIR spectrum is unaffected, while the X-ray 
flux decreases by a factor of 4, as expected from equations (1) and (2). 



4. Conclusions and Discussion 

The simultaneous detection of flares in the NIR and X-ray bands suggests that there is 
an intimate connection between these emission mechanisms. The fact that NIR flares have 
very steep spectra, whereas X-ray flares always have a much flatter spectrum, rules out the 
possibility of producing the NIR and X-ray photons together via synchrotron emission. The 
simplest model for production of X-ray is SSC model, but it faces several challenges. In this 
paper we have shown how the simultaneous observation of NIR and X-ray flares may be used 
to determine the source size, the magnetic field, and the distribution of electrons, should the 
NIR and X-ray emission be produced, respectively, via synchrotron and SSC. We have shown 
that the emission characteristics suggest an electron distribution cutting off at 7 cr ~ 100. 
The size of the source must be a fraction of a Schwarzschild radius in order to produce 
prominent X-ray flares. The rapid rise in the NIR emission suggests a magnetic field of a 
few tens of Gauss. A lower magnetic field is required should the source be Doppler-boosted, 
producing flares with soft NIR and hard X-ray spectra. 

We have demonstrated that the SA model we proposed previously is fully consistent with 
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the current observation of flares from Sagittarius A*, and with the transit-time damping ac- 
celeration taking the place of the parallel propagating waves we used before, the acceleration 
only depends on the turbulence energy density. However, depending on details of the flare 
excitation mechanism, the electron distribution below 7 cr can be quite different. This model 
therefore can be a powerful tool in probing the plasma energization and particle acceleration 
processes near the event horizon of the black hole. 

The rest frame acceleration time is slightly longer than the typically observed rise time. 
This suggests that a time-dependent treatment of the SA model is necessary unless the 
emitting plasma is Doppler-boosted by a boosting factor T ^> r ac /r r i S . However, our main 
conclusions still hold — current observations cannot yet distinguish the subtle difference. Nev- 
ertheless, such a calculation is clearly warranted with the acquisition of new data, particularly 
as flares continue to be observed simultaneously across the spectrum. 

It is also interesting to note that should the magnetic field be anchored to a slow large 
scale flow, it will exert a stress on the plasma. This stress causes transport of the angular 
momentum of the emitting plasma. In a Newtonian potential where the angular momentum 
density in a Keplerian orbit is given by L = m v n{GMr) 1 l' 2 , we have m p n(GM) 1 / 2 r _1 / 2 i' r ~ 
(VxB)-Br - B 2 r/R, or 




which is comparable to the Keplerian velocity at r = 3rs and is much higher than the 
radial velocity v r ~ 10 8 cm s" 1 suggested by the chirping behavior of the X-ray and NIR 
flares, indicating that the variation length scale of this external magnetic field should be 
much longer than R, or that most of the magnetic field is generated within the plasma itself. 
Should neither of these scenarios be viable, a large black hole spin would then be required. 

This model is also in line with our previous study of Sagittarius A*'s emission in the 
quiescent state (LPM04), and with our study of proton and electron acceleration on larger 
scales (~ 20rs; Liu, Melia, and Petrosian 2005). Clearly an appropriate modeling of the 
magnetic field structure in the accretion flow of Sagittarius A*, and in any possible outflow, 
is required in order for us to have a complete understanding of this intriguing object. With it, 
we should hope to test its predicted correlated variability over a broad range in frequencies. 

We thank Fred Baganoff for providing the most recent results of their multi-wavelength 
monitoring of Sagittarius A* and Huirong Yan for useful discussion. This research was 
partially supported by NSF grant ATM-0312344, NASA grants NAG5-12111, NAG5 11918- 
1 (at Stanford), and NSF grant AST-0402502 (at Arizona). FM is very grateful to the 
University of Melbourne for its support (through a Miegunyah Fellowship). 
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Fig. 1. — Left: Constraints on B and 7 cr of an SSC model for the NIR and X-ray flares in 
Sagittarius A*. The lower-left portion of the diagram is excluded because the SSC photons 
have a very soft spectrum there. The upper-right portion corresponds to models with very 
hard NIR spectra. For the region with small B, the acceleration time can be much longer 
than the rise time of the flare. The three steep dashed lines give the acceleration time 
(indicated), which is four times shorter than the synchrotron time at 7 cr . Viable models 
are therefore located in the central region bounded by the thick lines. The solid and long- 
dashed lines are for three different synchrotron and SSC spectra produced by electrons with 
an isotropic relativistic Maxwellian distribution. The corresponding spectral indexes a at 
1.4 x 10 14 Hz and 10 18 Hz (~ 4.1 keV) are indicated. For electron distributions different 
from the Maxwellian one, these lines will shift accordingly. The filled and the open circle 
are for the models indicated by the thick solid and dashed line in right panel of Figure 2, 
respectively, right: Same as the left panel, but for constraints on B, R, n, and / turb from 
the observed NIR (qnir) and X-ray (a x ) spectral indexes, the flux density at 1.4 x 10 14 Hz 
F N i R and 10 18 Hz F x . The source is approximated as a uniform sphere. 
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Fig. 2. — left: The steady-state electron distribution normalized at 4^y cr for values of 5 
indicated. The particles are injected at 7 in = 0.137 cr . The escape time is longer for smaller 
values of 5, and the spectrum converges to a curve very close to the one for 5 = 0.01 as 
8 — > 0. The solid line is for the Maxwellian distribution, right: A fit to the NIR and X-ray 
flare spectrum. The parameter values include a Lorentz factor 7 cr = 75 and a total number 
of electrons M = 3.8 x 10 42 . The emitting plasma is assumed to be a uniform sphere, with 
R and B as indicated. The data include the Chandra steady state flux (lower butterfly), 
and the peak X-ray spectral component of the 27 October 2000 flare (Baganoff et al. 2001). 
In the NIR, the data points represent peak fluxes during individual events, and the crosses 
(barely visible below the theoretical curves) are the NIR spectra measured during the flare 
of July 15th, 2004 (Eisenhauer et al. 2005). There are 4 theoretical curves: (thick solid) 
overall fit to the peak NIR emission of the July 15th event; (thin solid) the rising and decay 
phase spectrum of this NIR event; (dashed) a fit to the peak emission of the 27 October 2000 
Chandra flare; and (dotted) is the same as the thick solid curve, except that R = 0.8 r s . 
The data points to the left of the panel correspond to the peak flux of variable radio and 
mm emission (fluctuating on a time scale of days to weeks), that originates at larger radii 
(LPM04). 



